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Abstract 

In this paper we consider the problem of pricing a perpetual American 
put option in an exponential regime-switching Levy model. For the case 
of the (dense) class of phase-type jumps and finitely many regimes we 
derive an explicit expression for the value function. The solution of the 
corresponding first passage problem under a state-dependent level rests 
on a path transformation and a new matrix Wiener-Hopf factorization 
result for this class of processes. 
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1 Introduction 

Consider a riskless bond and a stock whose price processes {B t ,t > 0} and 
{St, t > 0} are given by 

Bt = exp (J r{Z s )ds\, S t = exp{X t ), X = x, (1) 

with r(-) > the instantaneous interest rate, Z a finite state Markov process 
and X = {X t ,t > 0} a regime-switching phase-type Levy process (that will be 
specified below in Section [2}. When X t is a Brownian motion with drift and 
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r(-) is constant, the model |T]) reduces to the classical Black-Scholes model (BS). 
It has been well documented in the literature that the BS model is not flexible 
enough to accurately replicate observed market call prices simultaneously across 
different strikes and maturities. 

To address some of the deficiencies of the BS model it was proposed to 
replace the geometric Brownian motion by an exponential Levy process, mod- 
elling sudden stock price movements by jumps. A substantial literature has 
been devoted to the study and application of Levy models in derivative pricing; 
popular models include the infinite jump activity models, such as the NIG [5J, 
CGMY [TO], KoBoL [5] and hyperbolic processes [15], and the finite activity, 
jump-diffusion models - see also Cont and Tankov [TT] for an overview. In the 
latter category, for instance, Kou [21] investigated the pricing of European and 
barrier options in the case of double-exponential jumps; Asmussen et al. [4 
considered perpetual American and Russian options under phase-type jumps. 

In a parallel line of research the BS model was extended by allowing its 
parameters fi, a and r to be modulated by a finite state Markov chain Z . The 
process Z models (perceived) changes in economic factors and their influence on 
the stock price. See Guo [T5], [16] for background on this regime-switching model 
and further references. In the context of option pricing, Guo [TH [J5] and Guo 
and Zhang [T7] obtained closed form solutions of European, perpetual American 
put and lookback options for a two-state regime switching Brownian motion; For 
the case of N states, Jobert and Rogers [TB] considered the perpetual American 
put and numerically solved the finite time American put problem. 

In the present study we consider the model (JTJ) which combines both the im- 
portant features of regime-switching and jumps, motivated by the observation 
that Levy models have been successfully calibrated to options with single, short 
time maturities whereas regime-switching models fit well longer dated options. 
The model (JTJ allows, at least in principle, for a flexible specification of the 
jump-distribution, since the phase-type distributions are dense in the class of 
all distributions on a half- line (see [U Prop. 1]). Under this model, we ob- 
tain explicit, analytically tractable results for the value function of a perpetual 
American put and corresponding optimal exercise strategy under this model. 
Guo & Zhang [TT] and Jobert & Rogers [T5] have shown that the optimal stop- 
ping time takes the form of the first-passage problem of X t under a level k(Z t ) 
that depends on the current regime Zt ■ We will show that the optimal stopping 
time still takes this form in our model and subsequently solve the corresponding 
first-passage problem. The solution of the latter rests on a path transformation 
and new matrix Wiener-Hopf factorization results, which extend the classical 
factorization results of London et al. [23] . The results also extend Asmussen et 
al. [1] who solved the first-passage problem across a constant level in the case 
of two regimes using methods different from ours. 

To value a finite maturity American put under the model ([T]) the solution 
of the perpetual American put problem may in principle be used as building 
block in an approximation procedure that we will briefly outline now. In the 
setting of the BS model, Carr [9] investigated the approximation of a finite 
maturity American put price by randomizing its maturity and showed, by nu- 



2 



merical experiments, fast convergence of this algorithm. The proposed maturity 
randomization resulted in an iterative evaluation of a series of related perpetual- 
type American options, a procedure which was extended to the setting of jump- 
diffusions by Levendorskii [55] . The idea is then to combine our solution of the 
first-passage problem with Carr's ideas to develop a pricing algorithm of finite 
maturity American put options under a regime-switching Levy model - we leave 
further exploration of this idea for future research. 

The rest of the paper is organized as follows. Section [5J is devoted to the 
problem formulation and the solution of the perpetual American put problem 
in terms of a first-passage problem. The solution to the first-passage problem 
under a state-dependent level is developed in Sections [3J— El Finally, in Section 
[3 the case of two regimes is considered in detail. Proofs that are not given in 
the text are deferred to the Appendix. 

2 Problem formulation 
2.1 Model 

Let the bond and risky asset price processes be given as in |T|) such that E[Si] < 
oo, where Z = {Zt; t > 0} is a continuous time irreducible Markov process with 
finite state space E° — {1, . . . , N} and intensity-matrix G, and X = {Xt, t > 0} 
is a regime-switching jump-diffusion given by 

X t =x+ f n{Z s )ds+ f a{Z s )dW s +Y, f l(z.=i)dJi(«)- ( 2 ) 
Jo Jo ieE° J ° 

Here 1b is the indicator of the set B, x E M, W = {Wt;t > 0} is a Wiener 
process, Jj = {Ji(t);t > 0} are independent compound Poisson processes with 
jumps arriving at rate Ai, and /i and a are real- valued functions on E with 
cr(-) > 0. The stochastic processes X and Z are defined on some filtered prob- 
ability space F,P), where F = {Tt}t>o denotes the completed filtration 
generated by (X,Z). The jump sizes of the compound Poisson processes Ji 
are assumed to be distributed according to double phase-type distributions, the 
definition of which we will specify below. We first briefly review the defini- 
tion of a phase-type distribution. A distribution F on (0, oo) is said to be of 
phase-type, if it is the distribution of the absorption time of a finite state Markov 
chain with one state d absorbing and the remaining states transient. One writes 
F ~ PH(a, T) if this Markov chain, restricted to the transient states, has gen- 
erator matrix T and initial distribution given by the (column) vector a. From 
Markov chain theory it follows that the density of F is given by 

f(x) = a'e Tx t, x>0, (3) 

where ' denotes transpose and t = (— T)l, with 1 a column vector of ones, is 
the vector of exit rates from a transient state to d. The class of phase-type 
distributions is dense (in the sense of convergence in distribution) in the class 
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of all probability distributions on (0, oo). Examples of phase-type distributions 
include hyper-exponential and Erlang distributions. See Neuts [24] and As- 
mussen [TJ [2J [3] for further background on phase-type distributions and their 
applications. 

An extension to distributions supported on R reads as follows: 

Definition 1 A continuous distribution HonWLis said to be of double phase- 
type with parameters p, a, T, (3, U , and one writes H ~ DPH(p, a, T, j3, U), if 
its density h is of the form 

h{x) = pf a ,T(x)l( x>0 ) + (1 - P)f/3,u{-x)l {x<0) , (4) 

where p S (0,1) and f a ,T, fp.u are PH(a,T) and PH(f3,U) densities respec- 
tively. 

For each i 6 E° the jump size distribution F{ ~ DHP(pi, af , T ; + ,a~ ,T~) 
of Ji is of double phase-type. Note that the positive and negative jumps of Ji 
arrive at rates := PiXi and \ l _ := (1 — Pi)\i and are distributed according 
to PH(a^ l ",T i + ) and PH(a^ ,Tf) distributions, respectively. 

The market with price processes (B, S) as specified above is arbitrage-free 
as there exists an equivalent martingale measure P*. Furthermore, there exists 
a F* that is structure-preserving (i.e. X is still of the form |2]| but with different 
parameters) - a proof of this result is given in the Appendix. From now on we 
will assume that X admits a representation j2J under a martingale measure P* , 
and we will write P for P* . 

2.2 Perpetual American put 

In the market ((!]) we consider a perpetual American put with strike K > 0, 
a contract that gives its holder the right to exercise it at any moment t and 
receive the payment K — St- From standard theory of pricing American style 
options in QjJ] it follows that, if Sq = s = e x and Z = i, an arbitrage-free 
price for this contract is given by 

V*(s,i)= sup E at , i [B- 1 {K-e x -)+], (5) 
rer 0jCO 

where u + = max{u,0}, 7o i00 denotes the set of F-mcasurable finite stopping 
times, and E x ,i[-] = E[-|Xo = x, Zq = i\. An optimal stopping time in ((5]) reads 
as 

T(k*)=mi{t>0:X t <k*{Z t )}, (6) 

for some function (or vector) fc* : E° — > M. This can be seen to be true as 
follows. In view of the fact that (S, Z) is a Markov process, the general theory 
of optimal stopping in Shiryaev |27j implies that an optimal stopping time r* 
in ([5]) is given by 

r* = inf{t > : V*(S t , Z t ) <{K- S t )+}. 



4 



That t* is of the form © is a consequence of the fact that V* (s, i) is a positive, 
convex and decreasing function that dominates (K — s) + with V*(0,i) = K. 
The latter follows in turn by the definition of V* and by observing that s i— > 
[K — sc Xt ~ x ") + is convex and decreasing, and that subsequently taking the 
expectation and the supremum over the stopping times r S To.oo preserves 
these properties. 

The next result presents the solution to the valuation problem of the per- 
petual American put in the market (fTj): 

Theorem 1 The value function in §5§ is given by V*(s,i) = V k *{s,i) where 

V k (s,i)=Kv 0>k (x,i)-v hk (x ) i), ieE°,s = e x , (7) 

with 

v b<k (x,i)=E X}i [B- 1 ik) e bX ^ 
An optimal stopping time in §E§ is given by ([6]) where k* — {k\, . . . , k* N ) satisfies 
lim^(e l ,j) = -e*J j £ E°. (8) 

2.3 First passage 

To solve the American put problem we will consider the first-passage problem 
of X t under the level k(Z t ), with k : E° — > K, which amounts to finding the 
function 

v b , k (x,i) = E Xj i [e- RT+bXT h (Z T )] , (9) 

where Rt = Jq a{Z s )As, a, ho : E° — > R + , and 

T = T(k) = inf {t > : X t < k(Z t )}. (10) 

Following [2111], the first step in the solution of © is to reformulate this problem 
as a first-hitting time problem for a related continuous Markov additive process 
A, called the fluid embedding of X. Informally, a path of A is constructed from 
a path of X by replacing the jumps of X by linear stretches - see Figure Q] An 
explicit construction is given in Section [3] 

A classical approach [20- 22] to solving the resulting first-hitting time prob- 
lem rests on a characterization of the laws of corresponding up- and down- 
crossing ladder processes - see Figure [5] London et al. [23] developed matrix 
Wicncr-Hopf factorization results for fluctuating additive processes (see also 
Rogers [26] for elegant martingale proofs and Brownian perturbations). By ex- 
tending the results of [53] [55] to our setting, we solve the matrix Wiener-Hopf 
factorization problem for the embedding A, in Section 3] 

To deal with the different ways in which first-passage in can occur (see 
Figure [3J, the Wiener-Hopf factorization is employed in Section [5] to calculate 
the distribution of the process A at the first moment of leaving a finite interval 
or a regime-switch, whichever occurs first. In Section [6] the solution to ([9]) is 
derived by combining the foregoing results. 
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Figure 1: Shown is a sample path of X until the first time T~ that X enters 
( — oo, k). The process A has no positive jumps and always hits a level at first- 
passage. 



3 Fluid Embedding 

Let Y be an irreducible continuous time Markov chain with finite state space 
ELS d, where d is an absorbing cemetery state, and denote by A = {At, t > 0} 
the stochastic process given by 



A t =A + s(Y s )dW t 



m(Y s )ds, 



(11) 



where s and m are functions from E U d to K with s(d) = m(d) — 0. The 
process A is the fluid-embedding of X if the generator of Y restricted to E is 
equal to Qq where, in block notation, 




Qa 



Here D a is an N x N diagonal matrix with (D a )a 
of appropriate sizes and, again in block-notation, 



(12) 

A; +aj, O are zero matrices 



±„±' 



Xfaf 



B ± = 



N N 



If 



G 



A 



t 




From the form of Q a it follows that E can be partitioned as E — E + U E° U E~ 
where E° is the state-space inherited from Z and E + and E~ are the states 
in which the path of A is linear with slope +1 or —1 which originate from 
the positive and negative jumps of X, respectively. Similarly, we will write 
Ei = Ef U {i} U E~ for the subset of E corresponding to the «th regime of X 
with corresponding embedded positive and negative jumps. The functions m(-) 
and s(-) are then specified as follows: 



. , [ (T, if j € E° f ■ 

• W = otherwise ifJe £ ° 



3.1 Path transformation 



l if j e 
i ifjeE 
1 if j e 



In Figure [T] it is illustrated how a path of A can be transformed to obtain a path 
of X. More formally, denoting by 

To(t)= f l(Y B eE«)ds and T ~ 1 (u) = inf{£ > : T (t) > u}, 
Jo 
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the time before t spent by Y in E° and its right-continuous inverse, respectively, 
it is not hard to verify that 

( A o T ~ 1 , Y o T 1 ) is in law equal to (X,Z). 

This implies in particular that the triplets (Tq(T), Af, Yf) and (T, Xt, Zt) have 
the same distribution, where 

f = f(k) = inf {i > : Y t e £° and A t < fc(Yt)} 

with fc : £ — ► R given by k(j) = fc(i) for j E Ei. It is not hard to verify that 
state-dependent discounting (or 'killing') at rate a(i) when Y t = i G E° can be 
included by replacing Q by the generator Q a for the vector a = (a(i),i £ E°). 
Thus, it holds that 



v b , k (x,i) =E x . l ]c hA fh(Y f ) l(f <c ) 



(13) 



with y now evolving according to the generator Q a and E^ x [-] = E[-|F = 
i,A n = x]. Here, 

C - inf{t > : y t ^ (14) 
with inf = oo, and h : E — > K is given by /i(j) = /io(«) for j 6 



4 Matrix Wiener-Hopf factorization 



The solution of the first-passage problem of the Markov process (A, Y) across 
a constant_leveJ^ is closely linked to the up-crossing and down-crossing ladder 
processes Y + , Y~ of (A, Y). These processes are defined as time changes of Y 
that are constructed such that Y is observed only when A is at its maximum 
and at its minimum respectively, that is, 



y (^ + ) 



where 



t+ = inf{s > : A s > t) 



and 



and 



Yr 



(15) 



T t = inf{s > : A s < t}. 



It is easily verified that the ladder processes Y + and Y~ are again Markov pro- 
cesses with state spaces E°liE + ^&nd E°UE~ , respectively. We will characterize 
the generators and Q~ of Y + and Y~ along with the initial distributions 
rj + , defined by 



0,i 



Y+ = j, r+ < C] for iEE-,jeE+iJE a (16) 
Y ~ = j, to < C] for i G E + , j e E~ U E° (17) 



Denote by Q(n) the set of irreducible n x n generator matrices (matrices with 
non-negative off-diagonal elements and non-positive row sums) and write V(n, m) 
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for the set ofnxra matrices whose rows are sub-probability vectors. Let V and 
E denote the \E\ x \E\ diagonal matrices given by diag(ra(i)) and diag(s(i)), 
respectively. The matrix Q a is called recurrent if its rows sum up to zero; other- 
wise it is called transient. By considering the process A at the subsequent times 
it visits a certain state, r £ E say, and noting that this defines a random walk, 
we have that in the recurrent case either sup t>0 A t — oo or lim^oo A t = — oo, 
Po,»-a.s. 

Write N,N + and N~ for the number of elements of E°,E + and E~ , re- 
spectively and let iV + = N + N+ and 7V " = N + N~ . We then define the 
Wiener- Hopf factorization of (A, Y) as follows: 

Definition 2 Let G+,C+, G~ and C~ be elements of the sets Q(N+), V{N~,N£), 
Q(Nq ) andV(N + ,Nq), respectively. A quadruple (C + ,G + ,C~ ,G~) is called 
a Wiener-Hopf factorization of (A, Y) if 

E(-G+,W+)=0 and 3 (G~ , W~) = O, (18) 

where, for matrices W with \E\ rows, 

Z(S,W) = ^ 2 WS 2 + VWS + Q a W (19) 

and W + and W~ are given in obvious block notation by 

II I () /, | and W I /,. (20) 





where Jo, I + and I are identity matrices of sizes N x N, N + x N + and 
N~ x N~ , respectively, and O denotes a zero matrix of the appropriate size. 

In the following result the Wiener-Hopf factorization of (A, Y) is identified: 

Theorem 2 (i) The quadruple , Q + , n~ , Q~) is a Wiener-Hopf factoriza- 
tion of (A,Y). 

(ii) The Wiener-Hopf factorization (rj + ,Q + ,r)~ ,Q~) is unique if Q is tran- 
sient or if Q is recurrent and A oscillates (that is, sup t A t — — inf t A t = oo). 

(Hi) IfQ is recurrent and lim^oo A t = — oo, there are precisely two Wiener- 
Hopf factorizations of (A,Y) given by ,Q + ,rf~ ,Q~) and 

(n+(I - In) + ,i, Q+(I - In), rf,Q~) , 

where /i is the left eigenvector of Q + corresponding to its largest eigenvalue, 
normalized such that \x\ = 1, where 1 denotes a column vector of ones. 



Proof: For I e R let <£>^ be given by the matrices 

&f(x) = W+ cxp(Q+(e-x)) and &J(x) = W~ cxp (Q~(x - £)) , (21) 



9 



where W + and W are given by (f2"0")) with C + = r/ + and with C = f] ■ The 
proof rests on the martingale property of 

M t + = U(Y tAT +,A t/ , T} -) and Mr=/-(Y tAT -,A tAT -) (22) 

with 

/+(«,a:) = e<$|(a:)/i+ and f_(i,x) = e<$7 (a;) (23) 

where /i+ and /i_ are ATg — and Nq — column vectors, respectively. Here, and 
in the sequel, ej denotes a (column) vector of appropriate size with element 
e i{j) = 1 if J = i and zero otherwise. To verify that M + is a martingale, 
observe first that, since Y + is a Markov process with generator Q + = and 
initial distribution rj + defined by (|16[) . Markov chain theory implies that for 
a; < i 

= 4$t(x)h. (24) 

The martingale property of M + then follows from (|23|) - (f24|) as a consequence 
of the Markov property of (A,Y). An application of Ito's lemma shows that 
/ = (f+(i, u),i 6 E) satisfies for u < t 

\s{iff"{i, u) + m(i)f'(i, u) + J2 Qij «) -/(«,«)) = 0, (25) 

where /' and /" denote the first and second derivatives of / with respect to u. 
By substituting the expressions (f2"Tj) - (f2"3")) into equation (|25p we find, since /i 
was arbitrary, that Q + and ry + satisfy the first set of equations of the system 
(fl~8l) . The proof for Q~ and 77" is analogous and omitted. The proofs of Theorem 
[2] (ii) , (iii) are deferred to the Appendix. □ 

Example (Gerber-Shiu penalty function) Let X t in |(5J) with x > be the 
surplus of an insurance company. Note that, for a = 0, N = 1, and in the 
absence of positive jumps, ([2]) reduces to the classical Cramer-Lundberg model. 
Denote by p = inf{f > : Xt < 0} the ruin time of X. Quantities of interest in 
this setting include the ruin probability P x ,i{p < 00) and the Gerber-Shiu [13j 
expected discounted penalty function which quantifies the severity of ruin by 
measuring the shortfall X p of X at the ruin time p. Under the model @ both 
these quantities can be expressed in terms of the functions defined in (|23|) . 
For instance, the probability of ruin in regime I G is given by 

P,,, (p < 00, Z p = £) = (x) ft, 

where fi(j) = 1 if j e Et and zero else and $q is given by (|2"Tj) with Q~ = Qq . 
More generally, for any non-negative function it on [0, 00) x E° the Gerber-Shiu 
expected discounted penalty function reads as 

E Xji [e- R »n(X p ,Z p )] =e<$o(x)ff, 
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where R p = f£ a(Z s )ds, <I> is given by (f2Tj) with Q = Q a , and g is the 
N~ -vector with elements 



'n(0,£) foriEE 
7r(— s, m)e' e e sTm t^ n ds for £ E E~ 



5 First exit from a finite interval 

The two-sided exit problem of A from the interval [k,£] for — oo < k < £ < +00 
is to find the distribution of the position of (A T , Y T ) at the first-exit time 

t = Tkj = inf{f > : A t i [kj]}. 

By considering appropriate linear combinations of the martingales M + and M~ 
defined in (|2"2")l we will now show that the two-sided exit problem can be solved 
explicitly in terms of (?7 + , Q + , Q~ ). To this end, introduce 

z+ = ( o^io ^ eQ+( ,_ fc)j z = ^ n~ Q ) e Q-<«-* )> 

and define the matrices 

* + (ar) = (w + e 0+( ' _a:) -W-e <r(a, - fc) 2 + ) (i-Z'Z+y 1 , (26) 

= (w-e 0- ^ - ^ -W + e Q+ V- x) Z-} (I-Z+Z')' 1 , (27) 

*°(s,x) = (e sx I -e si y+(x)J+ -e sk y-(x)J-)[-K( S )}- 1 , (28) 

where / is an identity matrix, are given in (|20p with C + = r/ + and C~ = rf~ , 
J 1 * 1 is the transpose of (f2"0"l) with C replaced by zero matrices, and 

K(s) = ^X 2 s 2 + Vs + Q a . (29) 

We will write e/^k e/^k 1 ^ we wish to clarify their dependence on k and I. 
The complete solution of the two-sided exit problem reads as follows: 

Proposition 1 Let h + , h~ and h) be functions that map E° U E + , E° U E~ 
and E to K. If 

s(i) 2 b 2 + m(i)b < -qa for all i E E, (30) 
then it holds for x E [k,£] and i EE that 

E xA [h+ (Y T ) l ( A T =*,r<C)] = e&+(x)h+, (31) 
E x>i [h- (Y T ) l ( ^ =fc , r<c) ] = e^-(x)h-, (32) 
E I ,«[e M c-/ l t(y c _)i ((<r) ] = e'^°(b,x)A h1 q a , (33) 

where q a = (—Q a )l, C * s defined by (|14[) and A^t is fie diagonal matrix with 
elements h'(j). 
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Proof: Define g+ and g_ by the right-hand sides of (|31[) and (|3"2"|) respectively. 
It is straightforward to verify from (|26| - (|27|) that it holds that 

h+(i) if x = £,i E E + U E° 
if x = k,i e E~ UE° ' 

if x = £,i e E+ UE° 
h~(i) ii x — k,i e E~ U E° 

In view of these boundary conditions and the fact that any linear combination 
of M + and M~, defined in ([22)1 . is a bounded martingale, Doob's optional 
stopping theorem gives that 

g+(i,x) = E x , l [g + (Y T ,A T )1 {T<C) ] 

= E„ [h+{Y T )l (AT=i;T<0 ] , 

where r = Tfe^. Similarly, it follows that 

g-(i,x) =E !E)i [h~{Y T )l {AT=k ^ T<c) ] . 

To prove the third identity, consider the map h* : E — > K given by 

h*(i)=E 0li [e'^-tf (y c _)]. 

By conditioning on the first jump epoch £ of Y it is straightforward to verify 
that 

~~ s(i) 2 s 2 /2 + m(i)s + gii ' 

where ^[fta] denotes the intensity of a transition i —> j [i — ► d]. After reorder- 
ing and writing the above expression in matrix form, it follows that 

K(s)h* = A ht Q„l, 

so that, for s satisfying s(i) 2 s 2 + m(i)s < —qu, i 6 E, 

h* = [-K(s)]- 1 A ht (-Q a )l. (34) 

In view of the strong Markov property and (|3"Tj) - Q320 , it follows that 

E xA [e sA <-l (c<T) ] = E xA [e sA <-]-E X!i [e sA <-l {(>T) ] 
= e sx h*(i)-E X:i [e sA -l (T<0 h*(Y T )] 
= e'i (e sx h* - e sk y-(x)J-h* - e se ^+(x)J+h*) . (35) 

Inserting ([34]) into ([35]) finishes the proof of ([33]) . □ 



g-(i,x) 
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(a) (b) (c) 

Figure 3: First-passage under the state-dependent level (fo, k\) takes place while 
Y G E° and can take place in two ways: A hits the level ki while Y — i, i = 1, 2 
(illustrated in (a) and (b)) or by a jump of Y (illustrated in (a), (b) and (c)). 
The different line styles of the paths of A correspond to the different states of 
Y. 



6 First-passage under state-dependent levels 

By combining the ingredients from the previous sections the first-passage func- 
tion Vk,b{x,i) of At under the level k(Yt) 

v blk (x,i) =E Xii [e bA ^h(Y f{k) ) l ( T (fe) <c)] , (36) 

with h : E — ► R+ and 

f(k) = inf jt > : Y t £ E° and A t < k(Y t )} 

can be explicitly expressed in terms of the matrix Wiener-Hopf factorization 
found in Theorem [5] For simplicity we will assume that the levels are ordered 
as k\ > &2 > • • • > ^jv (the general case of possibly equal levels follows by a 
similar reasoning). As the first-passage over k can only occur when Y is in E°, 
it follows that Y can cross the boundary k before A exits the interval [kj, kj—i] 
in two ways: either Y jumps into a state {1, . . . ,j — 1} or A hits the level kj 
while Y is in state j - see Figure [3] We are thus led to considering the processes 



where E J := E\ ufZi E, 



(with Y^ = Y). Clearly, the Y^ are themselves Markov processes with gen- 
erators Qtj) given by the corresponding restrictions of Q a ; in block notation the 
resulting partitions read as 
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for some matrix , where qk>) is the matrix of exit rates from the sub-space 
Ej . By the strong Markov property, Proposition [T] and Theorem [2l the value 
of Vb,k{x,i) can be expressed in terms of the unknowns Vb,k(kj,i). For these 
unknowns a system of equations can be derived by invoking smoothness and 
continuity properties of Vb,k above the barrier k. 

More specifically, denote, for some constants Cj(£), Dj(£), the vectors hj 

and ht ■ % -> K by 

hj{l) = 



respectively, and set 





iff G E°\{1, 


Dj-i(£) 


else 


e bk i h(j) 


if I = 3 


Cj(t) 


ill e E°\{1, 


DM 


else 



(j = 2,...,iV), (37) 



fet = (h(£),l € Ej). 



,N), (38) 



and 



We shall write , J , + , and , J , ° as shorthand for ^i,. fe *r k . k . ^ 
*(fc fc_i)' respectively, and denote by f{x+) and f(x—) the right- and left- 
limits of the function / at x. Then the following characterization of v holds 
true: 



Theorem 3 Assume that s(i) 2 b 2 /2 + m(i)b < —qu for all i € E. The function 
Vb,k is given by 



if x> ki, 



v b ,k{x,i) = < 



*+{x)h* + *j(x)hj + *°(6, as)A h tg^J , if J = 2, . . . , JV, 

kj < ~~~. x ^ — ^ ? 

(39) 

where $ = $ fei is given by (f2"3"| and /ij , /i, are specified by (j3"T|) - (J3HJ) wif/i 
Cj (£) and Z?j (f ) satisfying the following system of equations: 



(bl T: 



l£E°\{l,...,j}, (40) 

£e^\[^° ui; 7] . ( 41 ) 
_1 *7^(i), (42) 



where j = 1, . . . , N — 1, and ' denotes differentiation with respect to x. 



7 Example: two regimes 

To illustrate the results derived in previous sections we consider next the model 
([2]) in the case of two regimes. Suppose that Z is a Markov chain with state 
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space E° = {1,2} and transition matrix 



G = 



-qi qi 

92 -92 



and that X evolves as a Brownian motion with drift fi\t + <j{Wt when Z is in 
state 1 and as the jump-diffusion [i 2 t + a 2 Wt — Jt when Z is in state 2, with J 
a compound Poisson process with intensity rate A and exponential jumps with 
mean 1/a. Then the embedding of (X, Z) has state space E = {1,2,2*}, say, 
with corresponding transition matrix 



Q r = 92 -92 - r 2 - A A I . (43) 




We will consider the stopping time T{k\,k 2 ) f° r the three different configura- 
tions of the optimal levels: k\ < k 2 , k{ = k 2 and kf > k 2 . For x > max{fc 1; k 2 }, 
the value function of the put is determined by the generator matrix Q~ of the 
corresponding down-crossing ladder process, which we determine by invoking the 
matrix Wiener-Hopf factorization results from Section |H Noting that E + = 
and E~ = {2*}, it follows that Q~ = Q~ satisfies 

\y?{Q-) 2 + VQ- +Q r = 0, (44) 



where 




V 




Denoting by f3[&\ a right-eigenvector of Q r corresponding to eigenvalue 6 it 
follows by right-multiplying (|44|) with (3[9] that the matrix 

K{6) = i£ 2 2 + V9 + Q r 

is singular and K(8)(3[6] = 0. It is a matter of algebra to verify that 9 satisfies 
g(0) = where 

g(9) = F x (0) ((a + 9)F 2 (6) - X9) - qi q 2 (a + 9), (45) 

with 

F j (9) = l^e 2 +f, j e-q j -r j , j = 1,2. (46) 
The following result lists the properties of the roots of g(6) = 0: 



Lemma 1 Suppose that r\,r 2 > 0. Then g{9) has five different real roots which 
satisfy the distribution characteristics 9\ < 9 2 < 9% < < #4 < 9$. As a 
consequence, Q~ = Q~ has three distinct eigenvalues 9\, 9 2 and 9%. 



15 



Since the eigenvectors (3[9i\ corresponding to the different eigenvalues 6i, i — 
1, 2, 3, are linearly independent, the matrix Q~ explicitly reads as 

Q- = (9 1 P[9 1 ] 9 2 (3{9 2 ] 9 3 0[O a ]) {P[6 X \ 0[9 2 ] flfla])" 1 . (47) 

1. Case k = k\ = fc|. For x > k, the value function of the American 
put reads as 

/ K-e k 

W{e x ,i) = e' l e Q ~ {x ~ k) H{k) with H(k) = K - e k 



where Q is given in (|47|) and k solves 

-e fc = e'iQ-H(k) i = 1,2. 

2. Case k% > k{. To deal with the case that k\ < x < k^ and Z = 1, 
we note that, if the process Z is restricted to state 1, X is equal to a Brownian 
motion with drift, fi±t + <r\Wt, killed at rate q\ + r\. The generator matrices 
of this restriction of Z and the corresponding ladder processes, denoted by 
Q, — Q + , Q~, reduce in this case to scalars, given by 

Q = -qi - n 

and the positive and negative root of the equation 

^ofx 2 + fi!X - qi - rx = 0. 
The associated two-sided exit probabilities from the interval [k'l, k 2 ] are 

e Q-(x-k) _ e Q+(l-x) e Q- (£-k) 



*r(*) 



1 - c Q+(e-k) e Q-(i-k) 

e Q+{e-x) _ e Q-(x-k) e Q+(e-k) 
1 _ e Q'(£-k) e Q+(£-k) ' 



with k = k\ and I = fc|- Putting everything together shows that the value 
function of the American put in this case is given by 



W{e x ,i) 



where 



' e ^ e Q ^~K)H(k* 2 ), x>k%,i=l,2, 
H + {x) + Vi(x)[C - H+{k* 2 )] + %l(x)H-(kt), k{<x< k* 2 , i = 1, 

C 



H(k)=\ K~e k I, H-(k) = e k — 1 H+ (x) = K - c x 



1G 



using equation (|A.3[) . Here C is determined by 

W'{e k2 -,l) = e[Q-H{k 2 ) 
and the levels k\ and k 2 satisfy the smooth fit equations 

*?(k 1 )H-(k 1 ) + *t'fa)[C-H+(h i )] = -e fc \ 

e' 2 Q-H(k 2 ) = -e k \ 

where prime in the first equation denotes differentiation with respect to x. 

3. Case k\ > k 2 For k 2 < x < k\ and Z = 2, we are led to consider the 
Markov process Y^> with state space {2, 2*} and generator matrix 



Q(2 



-<?2 - r 2 - A A 



In this case it can be checked from (fl8|) . (|19|) . and ([20)1 with Q a replaced by 
Q(2), that — is a scalar given by the positive root of 

—a; + /i 2 .T + A — ■ — = q 2 + r 2 + A 
2 .t + a 



and that 



a/[- 



H2) 



We can calculate Q( 2 ) in a similar way as we calculated Q above. Writing Q + = 
Q%) an< ^ Q = Q(2) ^ ne corresponding matrices of two-sided exit probabilities 
from [k 2 , k^] read as 



'/ 



e Q+(£-a;) _ e Q-(x-fc) 



e Q+(«-fe) 



3 Q-(x-fe) _ e Q+(i?-x) Me Q-(£-fe) j _ e Q+{i-k) Me Q-(e-k) 



with k = k 2 and £ = k*, where 



M = 



1 - ( 1 )e c 

1 

rj+ 



"(<-*) 



The value function is 

e >. e Q-(x-K)H(k$), 



W(e x ,i) 



H+(x) + / *l(x)[C - H+(kt)] + ^ 2 (x)H-(k* 2 ) 



x ^ k ^ . 2 — 1 , 
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where / is the row vector / = ( 1 ) and 
'K - e k 



H(k)=\ C , H+(x) =K-e x ^ , H-(k)=e k 



D 



92 92 



a I i 



where we used again equation (|A.3[) . Here C and D are determined by the two 
linear equations 

W(e kl -,2*) = W(c kl +,2*), 

W'(e kl -,2) = W'(e kl +,2), 

and the levels k\ and k% satisfy the two equations 



^'(k 2 )H-(k 2 ) + *l'{k 2 ){C - H+ih)} 



- e k ^. 

92 

e'iQ-Hih) = -e kl . 
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Appendix 



A Proofs 

A.l Cramer Martingale Measure 

In this section and the next we present a construction of an equivalent martingale 
measure for the process (X, Z) and show how the parameters change under this 
change of measure. For a background of Markov additive processes we refer to 
Asmussen [3]. 

An important role in the construction of the change of measure is played by 
the process X a — {X a (t);t > 0} defined by 

X a (t) = / a(Z s )dX s 
Jo 

for some function a to be specified below. It is straightforward to verify that the 
process X a is still of the form @, but with changed parameters; its characteristic 
matrix is K a [s] = G + A a [s] , where A a [s] is the diagonal matrix with elements 
Ki(ciis). Write h and A for the Perron-Frobenius right-eigenvector and eigenvalue 
of -Ka[l], respectively, and define the candidate change of measure L = {L t ; t > 
0}by 

L{t) = e Xa ^- xt h{Z{t))/h{Z{Q)), (A.l) 
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with h(i) = hi the zth coordinate of h. If E[Sj.] < oo, then ftj(l) < oo for i € £7°, 
and a solution cii, % E E°, exists of the equations 

Kiffli + 1) = r(i) + Kiffli), (A.2) 

with Ki(ai) < oo, where «j(s) = logi?[e sX i] is the Laplace exponent of X\ — 
[lit + (7iWt + Ji(i). It is shown in the following result that the measure P* with 
Radon-Nikodym derivative 

dP* 

dP F = 

is indeed an equivalent martingale measure: 

Proposition 2 Suppose that E[Si] is finite. 

(i) The process L = {Lt,!Ft\t > 0} is a positive mean one martingale and 
P* is a probability measure; 

(ii) Under P*, exp (^Xt — J * r(Z s )ds^ = B^ 1 St is a martingale. 

In view of Proposition [2l the market JTJ) with price processes as specified 
in is arbitrage-free if E[Si] < oo. It is shown below in Proposition [3] that, 
under P*, X is still of the form Q but with changed parameters. The process 
B^ 1 St is martingale if the following restriction holds for the parameters of X 
(see also [H Sec. 2]): 

f^+ M (i)+A i (f ! i (l)-l)=r(i) ) (A.3) 

where 

Fi(s) = Pia'ii-sl - Ti)- l U + q^sl - U^u, 

denotes the moment-generating function of /j, the pdf of the jump-sizes of X 
in state i. 

Proof: of Proposition [D (i) Let (X, Z) be of the form ^ , with corresponding 
characteristic matrix K and suppose that g is a E°-vow vector. Asmussen and 
Kella [5] have shown that 

e bx *- ct g(Z t ) - e bx "g(Z ) - f e bx ^ cu g(Z u )du(K[b] - cl), (A.4) 

Jo 

is a row vector of martingales for c > and b such that the diagonal elements of 
K[b] are finite. Choosing in (|A.4|) the process X to be equal to X a , b = 1, c = A 
and g = h, it follows that e Xa ^~ xt h(Z t ) — hrZo) is a zero mean martingale. 
As /i is positive, the process L in (|A.1[) is thus a positive mean one martingale. 
The proof of (ii) can be found in the next section. □ 
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A. 2 Change of measure 

Proposition 3 Under P* , the process X is still of the form @ with a*(i) — 
a(i), 

= + ~ / 2/(1 ~ e^X^F^idy) 
Jo 

and with J* compound Poisson processes with changed jump rates 
A (+)* = A (+) F (+) [_ 7] and A (-)* = A C-> F<-> [7] 

and distributions of the positive and negative jumps of phase-type with represen- 
tations 

(a ( + ). j3; ( + ). ) = K ( +) a <))l; ( + ,a < ) ) flnd („(-)*,!;(-)*) = (at"-), 7^— )), 
respectively, where the parameters are transformed according to 

( a (+>7) )T (+,7)) = ( a (+)A + /F (+) [- 7 ],A; 1 r (+) A + +7/), (A.5) 
( a (-.7) )T (-,7)) = ( a (-)A_/F(-)[ 7 ],A: 1 T(-)A_ -7/), (A.6) 

where A + and A_ are i/ie diagonal matrices with respectively and 

on i/ie diagonal such that k + = (-7/ - T^)- 1 ^, fc_ = (7/ - T^)~H^ 

and I is an identity matrix of appropriate size. 

Proposition 4 Under P* 7 Z has intensity matrix G* with elements 
9tj = 9ijHj)/h(i), i ¥> h and 9u = -^2,9ijh{j)/h{i), 

where gij is the ijth element of G; 

Proof: of Propositions and [^J We first show how to find the characteristic 
matrix of X under P*. Denote by X a the process X a + bX, where b is such that 
the elements of the characteristic matrix of X a are finite, and let / be a function 
that maps E° to E. Applying Ito's lemma to e^ a(t) ~ At /i a (Z t )/(Z t ) shows that 

e^- xt h(Z t )f(Z t ) - e^h(Z )f(Z ) - f Xe^- Xu h(Z u )f(Z u )du 

Jo 

due ^a(«)-A« x 



/ 

■7 



(l{z.=i}^(i)/(i)'«i(o< + b) +Y,9iMi)f{3) ~ h(i)f(t))) 

is a P- martingale, where we wrote A = A a and h = h a . Since 

K(h a )i = (K a [l}h a )i = ^gi 3 (h a (j) - K(i)) + Ki(ai)h a (i), 
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it follows from taking expectations and rearranging terms that, in vector nota- 
tion, 

E a , l [L t e bX ( t h Zt } = U+ [ E 0ti [L u e bX ( u h Zu ](G* + A*[ft])du, (A.7) 
Jo 

where G* is as in the statement of the proposition and A* is the diagonal 
matrix with elements A*^ = Ki(ai + b) — m{ai). Writing F£[b] for the matrix 
with elements Eq i [e hXt lz t =j\ an d differentiating (|A.7|) with respect to t, we 
arrive at the matrix differential equation 

F*' [b] = F* [b] (G* + A*), F* [b]=I, 

where / denotes the time derivative. Solving this system shows that the charac- 
teristic matrix of X under P* is given by K* [b] = G* + A* . See pQ for a proof 
of (|A.5P and (|A.6[) ; the rest of the statements of (i) and (ii) directly follow from 
Proposition 3 in [25]. □ 
Proof: of Proposition\^u) From Proposition [3] it follows that the aj's have been 
chosen in such that, under P*, X % have cumulant-generating functions satisfying 
= Ki(a,i + 1) — Ki(cii) — ri, so that the characteristic matrix K* of (X — 
J Q r(Z s )ds, Z) under P* satisfies K*[l] = G* and 1 is an eigenvector of K*[l] 
with eigenvalue 0. Setting g = 1 in (|A.4|I and taking a = 0, it thus follows that 
the process e^ 4- -^ r ( z ^ ds is a martingale under P*. □ 

A. 3 Wiener-Hopf factorization 

Proof: of Theorem @ii) Now we turn to the proof of the uniqueness of the 
Wiener-Hopf factorization. To this end, let (Z + , G + , Z~, G - ) be a Wiener- 
Hopf factorization and define the function / as in (|2"Tj) , but replacing r/ + and 
Q + by Z + and G + respectively. Since (Z+,G + ) satisfies the first equation of 
(fT5|) . it follows from an application of Ito's lemma, that f(Y t ,A t ) is a local 
martingale that is bounded on {t < r^ + }, so that Doob's Optional Stopping 
Theorem implies that 

70»=E XJ [/(y tAT +,A 4AT/ )] 

= E x>j [f(Y+,A T +)l (T + < ^]+\im E Xij [f(Y t ,A t )l (T + =oo) ]. (A.8) 

By the definition of / and the absence of positive jumps of A, the first expecta- 
tion in (|A.8p is equal to f(j, x). Note that the second term in (|A.8[) is zero if Q 
is transient or Q is recurrent and sup t A t — +oo. Indeed, in the latter case, tf 
is finite a.s. whereas in the former case f X: i{Y t 6 E) converges to zero. Thus 
/ = / for all h and we deduce that G + = Q + and Z + = Similarly, one can 
show that G~ = Q~ and Z~ = rf~ and the uniqueness is proved. □ 
Proof: of Theorem [HJiii) Assume that Q is recurrent but A t — > — oo. As Q + 
inherits the irreducibility property of Q, it follows from the Perron- Frobenius 
theorem that the matrix Q + has a probability vector /jl as left-eigenvector with 
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its largest eigenvalue. Since the quadruple (r] + , Q + , r\~ , Q~) satisfies (fTSI) . it is 
straightforward to check that this remains the case if we replace (?7 + ,Q + ) by 
(r/ + (7 — In) + fx, Q + (I — We are left to show that these are the only two 

factorizations of (A,Y). As in the proof of Theorem [Hji), it follows that any 
factorization quadruple of (A,Y) must contain rf~ and Q~ . Letting (i] + ,G + ) 
and f(j,x) be as in the proof of Theorem^ we distinguish between the cases that 
G is recurrent or transient. In the latter case f{j,x) tends to zero if x — > — oo 
and we deduce from (|A.8[) that f — f and thus G = Q + and Z + = In 
the former case, we note that, as G + inherits the irreducibility property of Q, 
it has a unique invariant distribution v given by the left-eigenvector of G with 
eigenvalue 0. Thus f(j,x) converges to e'jlvh = vh as x — > oo. The right-hand 
side of (|A.8[) is thus equal to 

/0» = /0» +P XJ (r+ = oo)vh 

= f(j,x) + (1 - e^ + exp(Q+(> - x))l) vh. (A.9) 

By differentiation of (|A.9|) with respect to x, we deduce that G = Q + {I — If). 
In particular, it follows that v is a left-eigenvector of Q + . Since the Perron- 
Frobenius eigenvector is the unique eigenvector with the largest eigenvalue, it 
follows that = v and then also that Z + = — lfi) + /it, which completes 
the proof. □ 

A. 4 First-passage under state-dependent levels 

Proof: of Theorem\3\For brevity of notation we will drop the subscript and write 
v for Vb,k- Appealing to the strong Markov property, it follows that, for x > k±, 

v(x,i) = E x , i [v(ki,Y T -)l( T - <C )]=E Xli [hi(Y T -)l( T - < Q], 

where t~ = r^~, and for kj < x < 

v(x,i) = E x>i [u(fc i _i,y T )l (r<CiAx=fe; ._ l) ] 

+ E Kil [v{kj,Y T )l {T<(tAr=kj) ] + E x .i [u(A c _,y c _)l (c<r) ] . 

Invoking results from Proposition Q] yields that ([55)) is valid for some vectors 
hj , hj and hj . To finish the proof we have to show that the stated form of 
these vectors is correct. We start with noting that, by the structure of the 
process (A, Y), 

v{k h j) = e bk > and v(k Jt t) = e bk 'e' t {sl - Tryhj for I e EJ. 

Furthermore, we claim that i) is continuous. Indeed, from the Markov prop- 
erty it follows that for I £ m £ E° 

poo 

v(z,£)= / W (z + y,m)e^e T ™^ m dt/, (A.10) 
Jo 
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so that, in particular, it holds that v(-,£) is continuously differentiable on 
(k m ,oo). Similarly, it follows that v(-,l) G C 1 (fc m ,oo) for I G E+. The conti- 
nuity of v(-,£) for £ e E° follows directly from its definition. As a consequence 
it follows that the equation (|¥T|) holds true. Let £ > j,£ £ E° and consider 
v(x,£) for x G [fcj — e, fcj + e]. By a Feynman-Kac argument it follows that, on 
[kj — e, kj + (] for e > small enough (such that kj — e > fc^), £) is equal to 
the unique C 2 solution of the ODE 

^f" + mf - <Mf = 9, f{kj ± e) = v(k s ± e,*), 

for some continuous function g and some constant c(£). In particular, v{-,£) is 
continuously differentiable at kj and it follows that PO)) holds true. □ 



A. 5 American put 

Proof: of Theorem[l\The proof of this result follows a standard approach for solv- 
ing perpetual American option pricing problems. As argued above the optimal 
stopping time must be of the form ([6]). Therefore, the value function is given 
by Vk* for some vector k* € (— oo, \ogK) N . The vector k* can subsequently be 
found by optimisation. At this point we note that the condition (IA.3j) implies 
that for the embedding s(i) 2 /2 + m(i) < — qu is satisfied for all i € E, so that 
we can apply Theorem [3] Since, for fixed k t— * Vk(x, i) is continuously 

differentiable it follows that fc* satisfies 



9Vk , x 



= for all (x,i), j = 1,...,N. (A.ll) 



k=k* 



Consider next the finite difference [Vk(e kj+h ,j) — Vk(e kj ,j)]/h and note that it 
is equal to the sum 

Vk(& +h ,j)-V k+h (& +h ,j) V k+h (e k >+ h ,j)-V k (e k >,j) 

h + h • (A ' 12j 

Letting ft, | 0, it follows from (|A.11[) , that the first term converges to zero, while 
the second term converges to — e kj . Thus we see that the smooth fit equations 
((SJ hold true. By a martingale argument it also follows that Vk* = V* for any 
solution k* e (-00, log iT)^ of ©. □ 
Proof: of Lemma 1 Suppose first that 

a + \pi + y / Mi + 2(ri+ 9 i)cr 2 ]/(T 2 . 
From the definitions of g(0),Fi (6), and F2(0), we have that 

g(+oo) = +oo, g(-oo) = -oo, 5(0) = a[(gi + r 1 )(q 2 + r 2 ) - qiq-i] > 0. 
Note that Fi(0) has two different real roots #0,1 > > 6*0,2 with 



0,2 = -\ni + \h4 + 2(r + qi)o\\lal 



23 



we then have #0,2 7^ — ol. Also, 

0(0o,i ) = -5i52(a + #0,1) < 0, 
because qi,q2,a, 60,1 > 0. Therefore, we further have that 



9(00,2) = 



-qiq 2 (a + 60,2) < 0, if #0,2 > -a, 
-9132(0: + #0,2) > 0, if #0,2 < -a. 



I AaFi(-a; r) > 0, if #0,2 > -a, 
f J " |A<xFi(-a;r) < 0, if O , 2 < -a. 

In view of the intermediate value theorem the proof of the first assertion is 
complete. Since Q~ is a generator matrix, it is negative semi-definite and the 
final assertion follows. 
If 

2 1 o7I 1 ~ \^2l l„1 



a = [Mi + y/Hi + 2(n + gi)tr?]/<^, 

= —a is a root. By a similar reasoning applied to h{9) = g(9)/{9 + a) and 
h(— a) < it can be shown that h has four distinct roots (two positive and two 
negative ones). □ 
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